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Abstract 

> 

We discuss several problems concerning domain walls in the spin S Ising 
. model at zero temperature in a magnetic field, H/(2S), applied in the x di- 

ON , rection. Some results are also given for the planar (y-z) model in a transverse 



field. We treat the quantum problem in one dimension by perturbation the- 
ory at small H and numerically over a large range of H. We obtain the spin 



C . 

O ■ density profile by fixing the spins at opposite ends of the chain to have op 

O 

posite signs of S z . One dimension is special in that there the quantum width 
of the wall is proportional to the size L of the system. We also study the 
quantitative features of the 'particle' band which extends up to energies of 
order H above the ground state. Except for the planar limit, this particle 
band is well separated from excitations having energy J/S involving creation 
of more walls. At large S this particle band develops energy gaps and the low- 
est sub-band has tunnel splittings of order H2 1 ~ 2S . This scale of energy gives 
rise to anomalous scaling with respect to a) finite size, b) temperature, or c) 
random potentials. The intrinsic width of the domain wall and the pinning 
energy are also defined and calculated in certain limiting cases. The general 
conclusion is that quantum effects prevent the wall from being sharp and in 
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higher dimension would prevent sudden excursions in the configuration of the 
wall. 

PACS numbers: 75.30.Et, 71.70.Ej, 75.30.Gw, 75.10.N, 75.10.H 
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I. INTRODUCTION 



Recently, there has been a growing interest in the study of interfaces with non-trivial 
geometry. Such interfaces arise in a variety of situations including domain walls in random 
mag netsE§il§ fluid invasion in porous mediaEE, spreading on heterogeneous surfaces!, 
membranes and vesicles in biology,! and epitaxial growth in materials science.i'0 

In connection with such problems it is natural to ask whether quantum effects play a 
significant role. For static properties it is well established that in nonrandom systems, such 
as a spin 5* antiferromagnet, with only a nearest neighbor exchange interaction, J, there 
are various regimes. Near the critical temperature at Tn ~ JS 2 , thermal fluctuations are 
dominant. In the ordered phase, as long as Tj S < T < T c , quantum effects due to the 
finiteness of S are unimportant. For T < T c /S one is in a quantum regime, where the 
quantum statistics of spin waves and their interactions leads to dependences on 5* and T 
not present in the classical (S — > oo) limit. Heuristic arguments indicate that quantum 
effects could influence the nature of domain walls in spin systems. In the picture in which 
the trajectory of the domain wall is likened to the trajectory of a particle in space as a 
function of time, quantum effects cause a smearing out of the trajectory. For such a system 
analytic and numerical work is obviously very difficult. Accordingly, we have been led to 
carry out a program of analytic work for quantum domain walls in one dimension. For this 
purpose we consider domain walls at zero temperature T in the spin S Ising model (with 
nearest-neighbor exchange coupling J) in a transverse field, H, whose Hamiltonian is 

Hti = E S *W - ^-2 E S z (n)S z (n + 1) , (1) 

ZD n ZD n 

where S(n) is a quantum spin S operator at site n. 

To study domain walls in this model we introduce boundary conditions in which the 
spin at one end of the chain is fixed to be "up" and that at the other end is fixed to be 
"down." This model has some interest in its own right. For 5 = 1/2 its properties can, in 
principle, be related to those of the associated free Fermion model obtained via the Jordan- 
Wigne transformation. However, with domain-wall boundary conditions, this relation is 



not easy to implement. Accordingly, we approach the analysis of the properties of domain 
walls in this model via perturbation theory for H/ J <C 1 and more generally via numerical 
solution for the wave functions and energies of the ground state and the low-lying excited 
states. Needless to say, some of the properties of this model in one dimension can not be 
extrapolated to higher dimensional systems. However, in most cases, one can safely say 
which features can be so extrapolated and which can not. 

We may summarize our results for the one dimensional model as follows. The magne- 
tization profile has a width of order the length of the chain. The low-lying excited states 
comprise a manifold of "particle" states, which results when the center of the wall propagates 
from site to site. These results are easily understood within perturbation theory in H/ J. In 
the classical limit the width of the wall is of order the correlation length, i. e. it is of order 
a lattice constant, as long as one is far from the critical regime at H pa 2J, above which 
long range order disappears. In the classical limit and for small H/ J, we evaluate a barrier 
energy, analogous to the Peierls-Nabarro energy,0 which prevents the free motion of the 
domain wall. When the correlation length becomes very large this barrier energy becomes 
small and since it is harder to calculate in this limit, we did not attempt such calculations. 
It would be interesting to calculate this energy for a quantum system, but in the present 
case, since the width of the wall is of order L, we can say that this barrier energy vanishes 
in the thermodynamic limit. Also, in the large S limit, we find that the quantum hopping 
of the domain wall from one site to the next is actually analogous to a tunneling process, 
so that the hopping matrix element is not of order H, as it is for the low spin case, but is 
now of order H exp(— aS), where a = In 2. It would be of some interest to recover this result 
within a field theory of one space and one time dimension to describe such a quantum effect. 

Briefly, this paper is organized as follows. In Sec. II we define the class of models we will 
analyze in one dimension. In Sec. Ill we give numerical results for the magnetization profile 
for S — 1/2, S — 1 and S — 3/2 and compare these to analytic results we obtain using 
perturbation theory in the small H/J limit. In Sec. IV we give analytic results for large S. 
The classical results for S = oo are given both in the continuum limit (i.e., when the wall is 
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very broad) and in the discrete limit (when the wall is very narrow). Here we also analyze 
the quantum system for large S. In Sec. V we give various results concerning the nature 
of the energy spectrum in the presence of a domain wall. We give numerical and analytic 
results for the nature of the particle spectrum caused by the matrix element which allows 
the wall to hop from one site to the next. Finally, in Sec. VI we summarize the conclusions 
to be drawn from our work. 



II. DEFINITION OF ONE DIMENSIONAL MODELS 

The first one-dimensional model we consider is the spin S Ising model in a transverse 
field, with the Hamiltonian 

^ti = -^E^H"^E S z {n)S z {n + 1) . (2) 

ZD n=l ZD n=0 

This formulation of the model has the advantages that a) the domain wall energy (for 
H = 0) is independent of the value of S and is equal to J, and b) the mean-field transition 
temperature To (at which significant correlations begin to develop between neighboring spins) 
is of order J independent of S. To discuss domain walls, the spins S(0) and S(L+ 1) will be 
fixed by boundary conditions, as discussed below. Thus L is the number of "active" spins 
in the chain. We will often use the notation h = H / J. 

We will also consider the "yz" model in a transverse field, for which the Hamiltonian is 

«»« = - | E - ^ E ( s «( n ) s «( n + !) + ^("W n + !)) • ( 3 ) 

The most important difference between these two models is that the interaction term (pro- 
portional to e) in 7i yz allows the spins to tip away from the z-axis. For H = and with 
free-end boundary conditions, the "y-z" model with e = 1 has the continuous U(l) symme- 
try instead of the discrete Z 2 symmetry of the Ising model which results when e ^ 1. In 
fact, we recall the ground state phase diagramS of this more general model, shown in Fig. 1 
for spin S = |, where one sees (at zero temperature) the disordered phase (D), the ordered 
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ferromagnetic Ising phase (F), and the ordered oscillatory phase (O). In both ordered phases 
the spontaneous magnetization (S z (n)) is no n- vanishing (where (• • ■) denotes the thermody- 
namic average at temperature T). However, the (connected) two-point correlation function, 
G{R) = (S z {R)S z {0)) - (S z (R))(S t (0)), behaves as follows^ for R -> oo: 



G(R) ~ < 



R~ 1 ' 2 exp (-R/0 , phase D; 

R- 2 exp {-R/0 , phase F; (4) 

iT 2 exp [-2R/g) ■ ReBe iKR , phase O; 

where the correlation length £ is a known function of h and e in each of the three phases, 



B is a constant and cos if = y/i 2 /(4e). When going from phase F to phase O, there is a 
new diverging length A ~ 1/K, which is not simply related to a gap in the spectrum of Ti, yz 
(as is £), but indicates the wavelength at which the correlation function at large distances 
oscillates under an exponentially decaying envelope. 

To discuss the width and energy of a domain wall, we shall work with "up-down" bound- 
ary conditions, in which we require 

S Z (0) = -S Z (L + 1) = S . (5) 

We can then study the profile of the wall by evaluating 

M(n) = (0\S g {n)\0)/S , (6) 

where |0) denotes the ground state. As we discuss in more detail below, in order to obtain 
a spin profile for a quantum system, antiperiodic boundary conditions which introduce a 
boundary coupling (K/2)S Z (0)S Z (L + 1) cannot be used in a naive way. However, for a 
classical system such antiperiodic boundary conditions will prove convenient. The energy 
of the domain wall, E w , is defined to be E w = E — E p , where E is the energy with "up- 
down" or antiperiodic boundary conditions and E p that with "up-up" or periodic boundary 
conditions. 
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III. RESULTS FOR THE PROFILE FUNCTION FOR S = 1/2, 1 AND 3/2 



A. Numerical Results for the Ferromagnetic Phase 

We begin by presenting numerical results for the profile function M(n) for the ferromag- 
netic phase F. For general spin we are constrained to rather small systems, since the number 
of quantum states is (2S + 1) L for L sites (excluding the fixed boundary spins). So why not 
use translational symmetry to reduce the number of states? 

Indeed, that would be possible for the special case of antiperiodic boundary conditions 
(with K = J) mentioned in the previous section. These boundary conditions are peculiar 
since TC commutes with the modified translation operator T = a x (0)T where T is the usual 
translation operator and o~ x (0) changes the sign of <J Z (0). In addition, Ti can be decomposed 
(for any K) into a block diagonal form. The corresponding states are said to be in the even 
or odd sector, respectively, depending on their parity under spin reversal. Furthermore, 
M(n) = in that case and one should consider instead the quantity M(n) = (0\S z (n)\l) 
where |0) and |1) are the ground states in the even and odd sectors (with respect to the Ising 
symmetry) of the model, respectively. It is well-known that at least for periodic boundary 
conditions this definesEl a sensible finite-lattice approximation to the order parameter (see 
Ref. [H] and references therein). However, for K = J, even for L finite, there are two 
degenerate states |1), in the odd sector, one of which gives rise to the desired wall profile, 
while the profile of the other one is essentially flat. On the other hand, for K > J, one does 
get a wall profile, but its shape depends on the value of Kj J. We avoid these problems by 
choosing the "up-down" boundary conditions defined above in Eq. (||). Furthermore, since 
we are interested in situations far from criticality, the convergence of the finite-size data 
with L — > oo is usually quite rapid so that the moderate sizes achieved are sufficient for our 
purposes. 

The determination of the lowest eigenstates of the Hamiltonians 7Yti and 7i JZ using the 
Lanczos algorithm is fairly standardB. From the ground state |0) we calculate the local 
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magnetization m(r) 



m(r) = (0\S z (n)/S\0) , r = n/(L + 1) . 



(7) 



Our results pertaining to the ferromagnetic phase F are as follows. In Fig. 2, we display 



smallest L-values used, we first observe that the data collapse onto a single curve, which 
implies that the L considered are already sufficiently large as compared to the correlation 
length that they represent the L — > oo limit. Second, we note that the profile is quite wide 
and a continuous function of r. We find a similar data collapse for all the situations we are 
going to consider in the sequel. Next, we show in Fig. 3, that the overall shape of the profile 
is not a peculiarity of having spin S = \. Rather, we see that the profiles obtained with the 
parameters e = and h = 0.1 (Fig. 3A) and h = 0.5 (Fig. 3B) both spins S = 1 and S — § 
are very close to the one for S — \. Nevertheless, the finer details of the respective shapes 
become different as H increases. Finally, in Fig. 4, we show the effect of varying h in the 
S — \ case. Note that the profiles for a transverse field as small as h = 10~ 5 and as large as 
h — 0.1 are superimposed onto each other. The fact that the width w of the domain wall is 
proportional to the size of the system is a peculiarity of a one dimensional quantum model. 
In higher dimensions we would expect w to remain finite as L — > oo except possibly at a 
critical point where the correlation length diverges. 

We shall turn to a quantitative explanation of these findings below. Profiles for the 
oscillatory phase are discussed in subsection C. 



In order to gain some understanding of the results of the last section, we now present 
some simple approximate calculations. As it turns out, most of the physics of the problem 
for h small is conveniently described in terms of degenerate perturbation theory involving 
the manifold M. of the h — ground states. This procedure avoids the additional technical 



m{r) for h — 0.1 and e = and S 



| for various values of L. With the exception of the 



B. Simplified Calculations for the Small S Quantum Case 
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complications of the free fermion method for "up-down" boundary conditions. The boundary- 
conditions are that spin 5(0) is fixed to have S z = S and spin S(L + 1) has S z = —S. For 
S — 1/2 the manifold M. contains the states \n), for n = 1, 2, ... L + 1, where \n) denotes 
the quantum state (shown in Fig. 5) in which spins % with i < n have their ^-component 
of spin equal to +1/2 and those with i > n have their z-component of spin equal to —1/2. 
The only nonzero matrix elements of the Hamiltonian within the ground manifold are 

H n>n = E = ~(L-l) J, l<n<L + l; 

7^n,n+l = Wn+l,n = , 1 < U < L . (8) 

Considering only the manifold M., one finds the corresponding eigenvectors ip p , and eigenen- 
ergies E p (for p — 1, 2, . . . L + 1) to be 

L+l 

tP p = Cj2 Mnpir/(L + 2)}\n) , E p = E + Hcos[pir/(L + 2)] , (9) 

71=1 

where C is a normalization constant. Perturbative corrections to the energy will occur at 
order H 2 / J and to the wavefunction at order H/ J. We may calculate the profile function 
within the small h approximation: 

& sin 2 [A;W(L + 2)] - £Li sin 2 [A;pvr/(L + 2)] 



M p {n) = ^ p {S nz /S)^ p 



= 1 



E^isin 2 [W(^ + 2)] 

, ELiSin 2 [fcW(^ + 2)] 
£^sin 2 [A; W (L + 2)] 



i 2 (2n + l)sinf^-sin(^ 
(2L + 3)sin^-sin(^ • 

We set n = xL and work in the limit of infinite L. Thereby we find 

m p (x) = M p (xL) = 1 - 2x + — sin(2pa;7r) . (11) 

pre 

In Fig. 3, mi(r) is shown together with the numerically determined ground-state profiles 
m(r) for S — |, S — 1 and 5 = | and we find a nice qualitative agreement (even for 
moderately large values of h). 

The above calculation can be generalized to larger S. Consider a state in which the wall 
is between lattice sites. In this state let all spins to the left of the wall have S z = S and 

9 



those to the right have S z = —S. Note that it is possible to change the value of S z for 
either one (but not simultaneously both) of the spins adjacent to the wall without changing 
the unperturbed (h = 0) energy. For spin S — 1 since the nonzero matrix elements of the 
perturbation due to the transverse field are all the same, the problem is totally equivalent 
to a spin 1/2 chain of twice the length. This observation explains the fact that our results 
for S = 1/2 and S = 1 are indistinguishable. When S > 1, one has to account for the fact 
that to move the wall through one lattice constant involves matrix elements which depend 
on the initial and final values of S z . This case will be considered later. 

C. Profile Function for the Oscillatory Phase 

The discussion has been so far restricted to the ferromagnetic phase F. The oscillatory 
phase is distinguished from it by showing a non-monotonic decrease of the correlation func- 
tion [see Eq. (Q)]. How does this behavior manifest itself on the level of the magnetization 
profile m(r) ? 

To answer this question, we display in Fig. 6A, for S — |, m(r) for h = 0.01 and e = 0.5 
and for L even. Indeed, this profile is distinct form the ones observed previously in the phase 
F. First, finite-size effects are much larger than in the situations discussed before. Second, 
m(r) displays, at least for L finite, a step-like behavior and it looks as if the system was built 
out of hard objects each occupying two lattice sites. When h is increased, these composites 
soften until they completely melt at the transition towards the F phase. 

This picture is modified in interesting ways for L odd. In fact, the calculation of m(r) 
requires a little more care in this case. For L even, the system has a single well-defined 
ground state separated by a gap at least of order 0(L~ 3 ) (see Sec. V) from the excited 
states. That is not so for L odd. Rather, for h = but e ^ 0, one finds that the ground 
state is twofold degenerate. That is a new degeneracy which has nothing to do with the 
ferromagnetic ordering of the system. Even when h ^ 0, the two lowest states in Ai remain 
much closer to each other than with the other states, the latter one having gaps of order 
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0(L~ 2 ). This new (near) degeneracy implies that one must reconsider the calculation of 
the magnetization profile m(r). Rather than blindly using Eq. (|7]), we take the two ground 
states |0) and |0') to be nearly degenerate. Then we construct the matrix 



m r 



s- 1 



( (0\S z (n)\0) (V\S,(n)\0) ^ 



r = n/(L + l) (12) 



v (0|^(n)|0') (O'l^(n)IO') j 

and find its two eigenvalues m±(r). Each of those represents a magnetization profile and one 
of them is shown in Fig. 6B. The other profile is obtained by reflection around r = 0.5. The 
asymmetry in Fig. 6B, observed for L odd and finite, is easily understood in terms of the 
composite objects introduced above. Note that within the manifold A4, the y — y term in 
the Hamiltonian moves the wall through two lattice spaces. Since we are discussing h = 0, 
this is the only kinetic energy of the domain wall. Thus, if we write M. = M e + A4 Q , where 
A4 e (M. ) is the submanifold of state in which there are an even (odd) number of up spins, 
then there are no matrix elements between these two submanifolds. Furthermore, these two 
submanifolds are related to one another by spatial inversion (accompanied by S z — > —S z ). 
Thus, for L odd, the problem decomposes into two identical eigenvalue problems (hence 
the twofold degeneracy of the energies). The matrix to be diagonalized is exactly the one 
considered in Sec. III.B, but now for M — (L + l)/2 sites and with h — > —e. The two sets 
of eigenstates are thus 



2m - 1 > 



^ +) = ^E-n(|^)|2m>, (13) 



where C is a normalization constant. This eigenvector basis also renders the matrix ([12]) 
diagonal. Writing n = 2k — 1 and n = 2k, respectively, we find for the profile for the + set 
of eigenstates (indicated by a subscript " +" ) 

( (2k + 1) sin(27rp/(L + 3)) - sin((2A; + l)2np/(L + 3)) 



MI (2k - 1) = M+(2k) =1-2- 



(L + 2) sin(27rp/(L + 3)) - sin((L + 2)2tt P /(L + 3)) 

(14) 
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In particular, this reproduces the steps in mf(r) observed numerically. Now, taking the 
continuum limit, it is easy to see that the width of the terraces decreases as L — > oo and 
that fh(x) — > m(x) + 0(1/ L). We therefore suggest that the phenomenon observed in Fig. 
7 is a novel type of finite-size effect. As L — > oo, the step-like and asymmetric profiles found 
for L finite will converge towards a smooth and symmetric limit function. 

At least for h — > 0, this limit function appears to show a simple relationship with the 
profiles previously found in the F phase. Indeed, our data suggest the following. The profiles 
for L large as found for h = 0, e ^ are very close to those obtained for h = h e s(e), e = 0. 
For small e our analytic work shows this to be true with h e g(e) = e. For larger e we have 
the phenomenological result, 

Me) - v 7 ^ (15) 

Evidence for this is provided in Fig. 7, where some profiles found for h = and e 7^ (data 
points) are collapsed with profiles calculated with h = h e s(e). 

IV. PROFILES AT LARGE S 

In this section we obtain several analytic results in the large S limit. We first analyze 
the classical (S — > 00) system. Then we treat the quantum system with S large, but not 
infinite. 

To treat the classical (S = 00) system, we write 

S x (n)/S = cos0 n , S z (n)/S = sin0 n , (16) 

where is a continuous classical variable. Then the energy of the classical system is 

1 L 1 L - 1 

E = --H C °S 4>n - - J Sil1 0" Sin ■ (17) 

^ n=l ^ n=l 

For the classical calculation without a wall, we will adopt periodic boundary conditions, so 
that 0l+i = 0i and the last term in Eq. flT7| ) runs from n = 1 to n = L. In that case 
<Pi — 0eq = arccos h, where, for convenience, we set 
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H/(2J) = h/2 = h/h c = h , (18) 

where h c = 2 is the critical value of h above which long-range Ising order disappears. For 
the periodic chain the energy per spin, eo, is given by 

e Q = - l -J(l + ~h 2 ) . (19) 

In the classical calculations we will use two types of boundary conditions to generate a 
wall. The first type of boundary condition is used for the continuum calculation. Here the 
number of sites is infinite and hence we set fax = — oo) = — <f>(oo) = eq . This continuum 
calculation is valid when the angle <fi changes slowly with position, as it does for h « 1. 
The second type of boundary condition is that used in the discrete calculation. Here it is 
convenient to use antiperiodic boundary conditions in which spins at opposite ends of the 
chain are coupled antiferromagnetically, with K = J and we require 0_j = —fa. 

In any case, the ground state energy is found by minimizing E with respect to the 
fas. Apart from end effects this minimization yields the following conditions (for n = 
2, 3, 4, ... N — 1) which characterize the exact ground state, 

^(sin0 n+1 + sin0 n _i) = (H/2J)tan<f) n = htaxi(j) n . (20) 



A. Classical Continuum Limit 

For small h the angle will change abruptly at the domain wall and we will treat this 
case in the next subsection. As h increases (up to the critical value h = 1), n will become 
a smoother function of n. In that case, a continuum approximation for <f> n makes sense. To 
treat the continuum limit for h ~ 1, we set sin0 n = y(x) = sin<p(x), where x = na, whence 
Eq. (|20D becomes 

hy{x) = ^ \jl - y{x) 2 [y{x - a) + y(x + a)] . (21) 
The continuum limit of this equation is 
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hy{x) = Jl-y(x 



( \ , 1 2 d<2 y 



(22) 



The solution to this differential equation is reduced to quadratures: 



V2x ~ 1 (1 - hcos<p{x) + VI - h 2 sin<f)(x)\ 

'1 — ^2 \ cos — a J 



We now discuss the significance of this result. First we analyze the behavior for large 
x. At large \x\, 4>(x) approaches ±0 eq , where cos0 eq = h. Thus for large \x\ we set <j>{x) = 
±[0 eq — 6(x)]. Keeping only the dominant terms in Eq. (p3|) , we obtain 

eq (x = In M^- . 24 

This gives 

= 8 e- x 't , (25) 

where 

<5 = 2(1 - /?) exp (a arccos h/{V2£,)) (26) 
and £ is the correlation length, given by 

^^=- (27) 
a a/2 Vl - ft 2 

Thus, for h — > 1, 

e/a ~ ^(ft c - fc)"" , (28) 

with v = 1/2, h c = 1, and A 5 = 1/2. 

Next let us see the behavior of the solution near x = 0. To do that, differentiate the 
solution above with respect to x at x = 0: 

^=^+-L^, (29) 
a ax 1 — ft, ax 

where we used 0(0) = 0. This gives 

14 



V2 



(30) 



dx] x= o a 

To estimate the width, W, of the wall, we note that over a distance W, the angle <p varies 
from approximately — <p cq to approximately +<p eq . This reasoning indicates that 



dx. 



20 



x=0 



eq j 



(31) 



For h — > 1, we have cq ~ v^2 V 1 — ^, so that 



W « 2a(l - /i) 



_ M-V2 



4£ 



(32) 



Thus for h h c , W is of order £, as one might expect. 



B. Classical Wall At a Lattice Site 

In this and the next subsection we look at the discrete equations, which are the appro- 
priate ones for small h. Since we are dealing with a classical system, it is convenient to use 
antiperiodic boundary conditions in which spins at opposite ends of the chain are coupled 
antiferromagnetically, with K = J. 

First we consider the case when the center of the wall is at a lattice site. To treat 
this case, we consider a chain with an odd number of spins. Number the sites — N, —N + 
1, . . . , —1, 0,1, . . . ,N — 1,N and fix O — 0. Then 0_ n = — <p n and the energy E' for this 
chain of 2N + 1 sites is 

N-l i N 1 

E' = - J sin 4>n sin 0„+i - - J sin 2 <p N - H ^ cos n - -H . (33) 

n=l ^ n=l ^ 

The 4>n are determined by Eq. (^) for n = 2, 3, . . . N — 1, by 

- sin 2 — h tan 0i = , (34) 

and by 

1 1 

-sin07v-i + - sin0Ar — /itan07v = . (35) 
15 



To obtain the solutions for the 0's at small h we write 4> n — ir/2 — e n for e n <C 1. Then for 
I < j < N 

1-1 1 

/i cot ei = - cos e2 , /icotej = -[cosej-i + cose J+ i] , hcoteN = -[cosejv-i + cosejv] • 

Zi Zi Z 

(36) 

To obtain the energy at order h we only need to solve the equations correct to order h, in 
which case 

ei — 2h, e 2 = e 3 . . . = e N = h . (37) 
Then, to find the energy to order h, we may set sin0j = 1, and cos0j = (1 + Sj t i)h, so that 

-E'l J = (N - 1) + X - + h + 0(h 2 ) . (38) 
The wall energy E' w is given by the large N limit of 

E' w = E' - (2N + l)e , (39) 
where eo is the energy per site of the uniform chain. Thus 

E'JJ =-(N- 1 -)-h + (N+ 1 -) + 0(/?) = l-h + 0(h 2 ) ■ (40) 

C. Classical Wall Between Lattice Sites 

To treat the case when the center of the wall is midway between two lattice sites, we 
consider a chain of 2iV sites, numbered —N, —N + 1, . . . , —2, —1, 1,2, ... N — 1,N. We 
require that 0_ n = — <p n and we use the same antiperiodic boundary conditions as in the 
previous case. Then the energy is 

N-l i j N 

E" = sin s m 0n+i + - J sin 2 </>i - - J sin 2 <p N - H ^ cos <p n . (41) 

n=l ^ 1 n=\ 

As before, the wall energy, E'^ is obtained via 
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E' = E - 2Ne n . 



The (fi n are determined by Eq. (EH) for n — 2, 3, ... N — 1, by Eq. (p5|), and by 



-(sin 02 — sin (fix) = h tan (fix . 



(42) 



(43) 



To analyze these equations, set (fi n = eq — e n . Then, to first order in e n we have (for 
1 < n< N) 



/itan^eq - — 2 A e n = sin0 eq - ^ cos o ( , ,[(-„_, + r„ + . 



cos^ 



'eq 



or 



(44) 



(45) 



We expect an exponentially decaying solution for e n for n 1. So set 



-na/5 



(46) 



so that 



(47) 



from which we get 



/i 2 \ /i 4 



1 - h 4 



(48) 



For small h, we get £ « a/(2| ln/t|). For /t — > 1, we get £ = (a/2)(l — ft) -1 / 2 , in agreement 
with the continuum result. These results hold irrespective of the position of the center of 
the wall. 

We now analyze the solution near the wall for small h. To do that we first look at Eq. 
(|3|) for 0i when we set (fij = n/2 — ef 



h cot ei = - [cos e 2 — cos ei] . 



(49) 



But e 2 <C ex, so 
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So 



ei = (4/i) 1 / 3 . 



(51) 



From here on one should interpret h to mean \h\. Next, look at the equation for 02 = vr/2— t^- 

h cot t2 = - [cos €\ + cos e 3 ] ~ 1 . (52) 

In this way we find that e n = h, for n > 1. 

A more systematic approach shows that we can write the solution for the e's as 

Cl = (Ah) 1 / 3 F 1 [(Ah) 2 ^] , e n = 6 cq + A n (Ah) 2 ^ 7 ^F n [(A~h) 2 / 3 } , n > 1 , (53) 

where A n = 2^ 6 ~ 5n \ e eq = sin -1 h, and the functions F n are analytic functions such that 

F n (0) = 1. 

To calculate E", note that, up to order h 4 ^ in the energy, we may write 



sin 



h = cosd = 1 - ^/^[i^/?/ 3 )] 2 + l(4/ i ) 4 / 3 [F 1 (/ i 2 / 3 )] 4 . . . 



sin (fij = cos €j = 1 , j > 1 ; 

cos0! = sin ei = (4/i) 1/3 F 1 (/i 2/3 ) + 0(h) 

cos0j = sinej = 0(h) , j > 1 . 

Putting these evaluations into Eq. (f41~D, we obtain 

-£"/J = 1 - ^(Ah) 2 ' z [F t (h 2 l*)f + ^(4/ i ) 4 / 3 [F 1 (/i 2 / 3 )] 4 + (N - 2) 



(54) 



1 - (Ahf' z [F x (h 2 ' 3 )f + \Ah)^F l (h 2 ' z )Y\ + - + 2/ i [(4/i) 1 / 3 F 1 (/ i 2 / 3 )] 

3 . 2 



= iV-l + (3/8)(4/i) 4/3 + 0(/i 2 ) . 
where we used the definition that -Fi(O) = 1. Then the wall energy is 

E'JJ = — (N - 1) - (3/8)(4/i) 4 / 3 + N = 1 - (3/8)(4/i) 4 / 3 . 



(55) 



(56) 



Note that this result is not identical to the case when the center of the wall is at a lattice 
site. The physics of this result will be discussed in more detail in the next section. 



D. Quantum Domain Wall for Large S 



We now consider the quantum chain for large S. For a chain of L spins of magnitude S, 
there are 2LS + 1 states which have the same h = ground state energy and which therefore 
must be treated within degenerate perturbation theory. Numerical studies show that for 
S > 1, this manifold of states splits into bands which are slightly different depending on 
whether 2S is even or odd. For S half integer there are 2S bands, each having L states 
except for the middle band which has L + 1 states. Approximately, the centers of these 
bands are at energies 

\(L-l)J+±mH, (57) 

where m assumes the values 0, ±3/2, ±5/2, . . . ±5. For integer S one has 2S — 1 bands, 
each having L states except for the middle band which has 2L + 1 states. In this case the 
centers of the bands are given by Eq. (j57|) , but m assumes the values 0, ±2, ±3, . . ., ±S. 
As an example, in Fig. 8 we show the density of states for S = 3 within the manifold M.. 
As Fig. 8 illustrates, the bands become wider and the gaps between bands become smaller 
as one approaches zero energy (about which the levels occur symmetrically). A qualitative 
explanation of the occurrence of such band gaps is as follows. If the spin 1/2 case is likened 
to a hopping model, then the spin S case is analogous to a hopping model in which there 
is a periodically variable hopping matrix element. This problem is therefore analogous to 
that of an electron in a weak periodic potential, where one knowJll that even for arbitrarily 
weak potentials one has band gaps. 

We will now demonstrate the existence of the band gaps and obtain quantitative in- 
formation on the band widths for ft C 1 as follows. We will give a construction for the 
outermost (and narrowest) sub-band, for which m = S in Eq. fl5"7|). Consider the subman- 
ifold of states, M.s, which contains the states, shown in Fig. 5, |p) (for p = 1,2, .. . ,L). 
These states are defined so that all spins (if any) having i < p have Si Z = S, all spins (if 
any) having i > p have S iz = —S, and S px = S. One may verify that the state |n) is in the 
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subspace M. of states of minimum energy eigenstates, Eq, of the h = Hamiltonian, i.e., 
that J/(2S 2 ) Yli=o S iz S i+ljZ \n) = E \n). The states \n) are not orthogonal to one another, 
but do have an overlap which is small for large 5* and is given by 

(p\p+l) = (S z = S\S x = S) 2 = T , (58) 

for < p < L. RoseS shows that 

(S, = S\S X = S) = 45(0, tt/2, 0) = 2- s , (59) 

where D is a rotation matrix. Therefore, we use states which are orthonormalized to leading 
order in r, 

|n) = |n) - (r/2)\n + 1) - (r/2)\n - 1) , (60) 

where we interpret |0) and \L + 1) to be zero. 

We now work only to first order in r. Then one notes that unless n — l<k<n + l, 
hSkx\ n ) has zero overlap with the subspace M.. To see this, note that (S z = m\S x \S z = 
m) = 0. Therefore we write, correct to leading order in r, 

(n | H\n + l) 

= (n\H\n + 1) - (T/2){n\H\n) - (r/2)(n + l\H\n + 1) 

= ~2S ^ ^n|5 te |n + l)-(r/2)H5te|n)-(r/2)(n + l|5 te |n+l)J 

= -t^HS^. + S'n+i^ln + 1) - (r/2)(n|S , nx |n) - (r/2)(n + 1| S^^n + 1)^J 

= (^25(n|n + 1) - (Sr/2)(n\n) - (Sr/2)(n + l\n + 1)^ 

= --Ht = -H2- 23 - 1 , (61) 
2 



where H = H-E . Also 

(n|W|n) = (n|W|n) + 0(r 2 ) = -^(n^n) = . (62) 

Thus we expect this subspace of states with S x = S to form a sub-band centered at energy 
£?o — if/2 with a width determined by the hopping matrix elementS t(S) = — H2~ 2S ~ 1 . 
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Thus the width of this band should be 4|t| = H2 1 ~ 2S . Numerically, for S = 5 we found a 
band width of 0.002if in excellent agreement with our calculation. We thus have an analytic 
description of the outermost two sub-bands corresponding to S x = ±S. We did not consider 
an analysis of the remaining inner sub-bands. Note that the results shown in Fig. 8 are 
not for asymptotically large S. We did check that when S is large enough, the outermost 
subband does become symmetric, as one would expect for a one-dimensional density of states 
corresponding to Eqs. (|6l|) and fl62]). 

In a sense, one sees from this calculation that for the wall to move one lattice spacing, 
it must tunnel through the phase space of spin states indicated by the result of Eq. (p9|). 
In this connection a nice analogy has been suggested by Stinchcombe.il! Consider a rotor 
in a strong cos (29) potential. On knows from exact solutions of the associated eigenvalue 
problem that one has a limit in which one has harmonic oscillator levels in the bottom of 
the potential well, all of which are doublets with a small tunnel splitting.0 It is tempting 
to think that one can provide a field theory in one space and time dimension in which the 
single potential minimum of the rotor becomes a periodic potential in a lattice. Then the 
bands would reflect the tunnel splitting. 



V. NATURE OF THE SPECTRUM 

A. L~ 2 Spectrum in the Ferromagnetic Phase for Small S 

So far, we have considered the form of the profiles m(r) and have seen that at least 
qualitatively, their form can be understood from simple degenerate perturbation theory. 
The manifold M. considered consists of those states which give the lowest energies in the 
h — > limit. For the case of S = ~, for h = these are the states shown in Fig. 5, namely, 

ITT---U) It - • - U • - • i> IU---U) (63) 

and for L sites, there are L + 1 of them. For S = 1 one starts at h = from the states 
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|TT---TO|) |t---tO|...|) ITOI---U) (64) 

and one gets a subspace of (L + 1) + L = 2L + 1 states. All these states have in common 

that they contain a single wall. States which do not belong to this "one-particle subspace" 

M. have a large gap with the ground state. It is now remarkable that this clear separation 

of the subspace M. from all other states remains intact even for finite values of h. The 

subspace M. has a simple analogy for (anti)periodic boundary conditions. In these cases, Tt 

commutes with the (modified) translation operator T (or T) and is thus decomposed into 

block diagonal form, the blocks being labeled by the eigenvalues of T (or T). M. corresponds 

to the set of lowest eigenstates of the blocks of Tt. 

While the gap of the lowest excited state outside A4 with the ground state is finite (e.g. 

1 — h for an d when e = 0), the gaps inside M. scale with L~ 2 . If g(i) = Ei — Eq 

denotes the i th gap, we expect the scaling g(i) ~ L~ 9 with 9 = 2, as obtained in Eq. ([|). 

Finite-size estimates for 9 are then obtained from, with N = L + 2 

]n(g(i) N+1 /g(i) N ) 
° L = HN/(N + 1)) (65) 

The extrapolation towards the limit iV — > oo was carried out with the BST extrapolation 

algorithmB. We illustrate the convergence of the finite-size data with a few examples for 

S = 1 in table [IJ The convergence towards 9 = 2 is even clearer for S — |. 

We now consider the finite-size amplitudes 

a(i) = N 2 g{i) . (66) 

The motivation for this comes from systems being precisely at a critical point. In that 
case, it is known that for one-dimensional quantum chains finite-size amplitudes defined as 
Ng(i) are related^ to the anomalous dimensions of the scaling operators of the model and 



this leads to simple patterns of the spectrum of the amplitudes (see e.g. Ref. for more 
information). Although we are not working here at a critical point, we observe a simple 
structure of the excited states within Ai for L — > oo out of the critical region. Namely 
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r(z - 1) = o(z - l)/a(l) = i(j - l)(i + 1) (67) 

Evidence for this is presented in table || for S = | and table |T| for 5 = 1 which give the 
BST extrapolated^ estimates for the a(i). It is remarkable that the same expression should 
hold for the gaps ratios of both the S = | and the 5 = 1 models, even though h is not very 
small and/or 6 ^ 0. For h small and e = 0, this result can be recovered from the energies 
(§) in lowest order of perturbation theory. 

B. L~ e Spectrum in the Oscillatory Phase for Small S 

Turning to the oscillatory phase, we note that the spectrum of the states within the 
subspace M. is different from that in the F phase. Here there are important differences 
between L even and L odd. 

For simplicity, we restrict attention to spin S = \ and first concentrate on the line h = 0. 
The structure of the levels we are going to find is shown in Fig. 9. For L even, the spectrum 
of the low-lying states corresponding to M. is grouped into doublets. While the doublet 
splitting is of order 0(L~ 3 ), the gaps between doublets are of order (D(L~ 2 ). For L odd, on 
the other hand, the energy levels are doubly degenerate, with gaps between them of order 
0(L- 2 ). 

To see this, we look for a scaling behavior of the energy gaps g(i) = Ei — Eq ~ L~ e and 
define finite-size estimates for the exponent 9 using (|S5|) . In table [V], we give our results for 
e = 0.02. For L even, estimates of 9 from the lowest six gaps (i — 1, . . . , 6) are given. (For 
some of the higher gaps, the finite-size data for L small are not yet in the asymptotic regime 
and are thus discarded.) We observe that the lowest gap scales with an exponent 9 = 3, 
while all the higher gaps scale with 9 = 2. For L odd, on the other hand, each energy level 

is doubly degenerate even for finite L, i.e. g(2j) = g(2j — 1), j = 1, 2, Estimates of 9 for 

i = 2, 4, 6 are given in table [V] and we find 9 = 2 throughout. Similar results are also found 
for e = 0.1 and 0.5. 
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Next, we study the finite-size amplitudes a(i) = N 2 (Ei — Eo), with N = L + 2. Estimates 
for the lowest six amplitudes extrapolated to L —>■ oo are shown in table [VTj , obtained for L 
even or odd, respectively. We find that for L even, the estimates for pairs of amplitudes are 
quite close to each other and are consistent with being equal. Furthermore, their numerical 
values are near to the ones found for L odd. Taken together, the present data suggest 
the simple picture that in the L — > oo limit, the amplitudes a(i) should become doubly 
degenerate and independent of the evenness or oddness of L. Finally, when looking for the 
amplitude ratios a(i)/a(2), we find the following pattern 

0,1,1,- -5,... (68) 



which is consistent with the very same structure (67) of the low-lying amplitudes found on 
the e = line. This finding is certainly consistent with the relationship between the order 
parameter profiles observed earlier (see Fig. 7) between systems on the h = and e = 
lines. 

These findings can be reproduced from the perturbative arguments of section III.B when 
e is small. For L odd, the eigenvalue problem then reduces to one treated there, but with 
[L + l)/2 sites and h — > — e and each state being twofold degenerate. For L even, we get 
two distinct problems, one with L/2 sites and one with L/2 + 1 sites. Using the expression 
(§) for the energies, we find pairs of levels with a splitting of order 0(L~ 3 ) between them. 

Qualitatively, the same structure persists throughout the O phase. We illustrate this 
in Fig. 10, where we show for e = 0.25 and L = 10 sites the dependence of the first few 
eigenvalues on h. Because L is not very large, the statements we shall make about the 
spectrum apply most strongly to the lowest levels and less well to the higher levels. We see 
that inside the O phase, the lowest energies occur in pairs which oscillate around another. 
In the ground state, there are for L sites (L + l)/2 level crossings for L odd and L/2 
level crossings for L even, respectively. These level crossing do not continue into the F 
phase. Finally, we point out that at the F/D transition, the energies combine into a f|) 
representation of the c—\ Virasoro algebra, as predicted from conformal invariance0i3. In 
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agreement with that prediction, we observe in Fig. 10 close to the critical point h c = 1 + e an 
(approximately) equidistant level spacing with the degeneracies 1,1,1,1,2,.... (Actually, 
because of finite size effects, this structure is realized at h — 1.15 rather than at h — h c — 
1.25.) 

This kind of level crossings in the ground state of finite-lattice quantum system has 
been first investigated for periodic boundary conditions^!. There is was noted that the level 
crossings always occurred between the Z2 even and odd sectors of the model. Here we 
find that the level crossing persist even if the Z2 symmetry is broken by our chosen "up- 
down" boundary conditions (|5]). Applying finite-size scaling!! to the location hk of the level 
crossings, we expect (but did not test) a scaling hk(L) — h* ~ L~ x l v . For periodic boundary 
conditions it is known!! that v = 1/2. [h* = \f\k gives the O/F transition line). 



C. Intrinsic Wall Width 

It is obvious that the fact that the width of the wall for the quantum model is of order 
the system size is a result of a quantum superposition of states each of which have a narrow 
wall. It is tempting, therefore, to introduce a measure of the "intrinsic" width of the wall, 
which is the minimum width obtainable within the subspace of states under consideration. 
We emphasize that this concept depends on the subspace of states begin considered. As 
we vary e, the strength of the y-y coupling, it seems plausible that the character of the low 
lying "particle" states may vary. Thus, it would be of interest to introduce a measure of the 
intrinsic width applicable to this case. For this purpose consider the quantity Q defined as 

Q = (2S)- 1 M 8 ,x neM (n\S z (t) - S z (i + l)|n) , (69) 

where the maximum is to be taken over all possible quantum states \n) in the subspace Ai. 
(As long as the site i is not near an end, Q will depend only weakly on i.) If the states consist 
of sharp walls, then Q will be unity. So we use this quantity to define the intrinsic width, 
Wi, via Wi ~ 1/Q- This type of definition is somewhat similar to the inverse participation 
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ratio introduced to characterize localized and extended states of an electron in a random 
potential.!! Here we take the subspace M. to be the "particle" subspace of states. This 
definition is appropriate as long as the subspace of particle states can be unambiguously 
separated from the continuum. This requirement suggests that h <C J. But, we can also 
use this idea when e is nonzero. 

We should mention that the maximization required in Eq. (p9|) is easy to carry out. One 
simply forms the matrix, P, (in any representation) where 

P n , m = (2S)- 1 (n\S z (t) - S z {i + l)|m) , (70) 

and Q is the largest eigenvalue of P. When Q is less than unity, one also obtains information 
on the shape of the wall via the eigenvalue spectrum of P. For instance, if the wall has width 
3, and if the three largest eigenvalues of P are degenerate, then one would conclude that 
the intrinsic wall profile is linear. 

One context in which this concept provides some information is when we consider the 
large S limit of the transverse Ising model. In that case we have already seen that the low 
lying particle spectrum splits into bands corresponding roughly to values of S x . When we 
apply Eq. ( ]6"5| ) taking the subspace M. to be the lowest band of L states with S x = S, we 
expect to find Q = 1/2. We in fact verified this expectation in our numerical calculations. 
They showed that one obtains two nearly degenerate maximal eigenvalues of magnitude 
essentially equal to 1/2. There are two eigenvalues because there are two wavefunctions 
which can maximize Q since the minimum intrinsic wall width is two lattice spacings. We 
see then, that the prescription of Eq. (^9|) coincides with the perturbative calculations for 
small h. 



D. Barrier to Domain Wall Motion 

First we summarize the result for the classical (5* — > oo) system. Note that E' w < E'^. 
for small H . So to move the wall through one or more lattice constants requires an energy 
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E" — E' = J 



(71) 



We have numerically studied the energies of the classical domain wall when its center is 
either at a lattice site or midway between lattice sites. We find that the barrier energy A 
does not change sign as H is increased as might be suggested by Eq. (|7T|), but rather A 
decreases monotonically to zero as H is increased towards its critical value of unity 

We carried out similar calculations for quantum systems with S — \ and S = 1. We 
found that the difference in energy when the center of the wall moved from a lattice site 
to a point midway between lattice sites was too small to be accessible to double precision 
arithmetic. Since, for a chain of length L < 20, we can not imagine a parameter on such 
a scale, we assert that the quantum barrier energy is zero. Practically, we performed this 
test as follows. For sufficiently large lattice sizes L, we expect for the ground state energies 
Eq(2N) (i.e. wall between two sites) and Eq(2N + 1) (i.e. wall at a site) 

E (2N) = A-2N + B + ... 

(72) 

E (2N + 1) = A ■ {2N + 1) + B' + . . . 

where A is the ground state energy per spin and independent of the boundary conditions. For 
£ = ~, A = — ~(h+ l)E(\^Ah/ (h + 1)) where E is a complete elliptic integral^. The desired 
wall energy difference should be proportional to B — B'. This in turn can be estimated from 
the quantities 

p N = E (2N + 1) + E (2N - 1) - 2E (2N) ~ 2(B' - B) + . . . 

q N = E (2N + 2) + E (2N) - 2E (2N + 1) ~ 2(B - B') + . . . (73) 



and we note that pn and qN should have opposite signs. We give in table [TV] estimates 



for both pjy and for both S = | and S = 1 for h = 0.1. We chose this value of h to 
be small enough that is was for from criticality but large enough so that \B — B'\ should 
be easily observable if nonzero. Taking the entries at face value, the wall energies must be 
less than 0(1O -5 ). Moreover, since p and q should have opposite signs, the values given 
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in table (TVl appear to be entirely due to the leading finite-size corrections and we find that 



Pn, qjq decrease rapidly with increasing N. So we get an upper bound 

barrier energy ~ \B - B'\ < C(10 -8 ) (74) 

for both S = | and 5=1. This result is not surprising. In the limit when the width of the 
domain wall diverges (as it does for L — > oo), the barrier energy should vanish. Whether it 
vanishes algebraically or exponentially would require further numerical or analytic analysis. 



E. Crossover From Quantum to Classical Behavior 

Here we discuss the various regimes which exist at at large S. It is clear that bands of 
energy levels exist on various scales. It is useful to discuss the consequences of these energy 
scales in terms of scaling functions. For notational simplicity, in this subsection F(x) will 
denote a scaling function of x which, in general, will be different in different appearances. 

We first discuss the thermodynamic properties of the wall. For this purpose we will 
consider the entropy associated with the wall, S w , defined as the entropy with "up-down" 
boundary conditions minus that with "up-up" boundary conditions. We emphasize that the 
scaling of the bulk thermodynamics, i.e., the extensive part of the thermodynamic potentials 
is independent of the boundary conditions and only depends on the bulk thermodynamic 
variables.il At temperatures small compared to J/S the entropy with "up-up" boundary 
conditions will be zero, because the ground state is nondegenerate. The entropy with "up- 
down" boundary conditions will be that due to the band of states with just one domain wall. 
First there is an extreme low-temperature limit in which the entropy is vanishingly small: 

S w ss ; T < H2 1 - 2S L~ 2 , (75) 

where we always set ks, the Boltzmann constant equal to unity. In this limit only the 
lowest single state of the tunneling band found in Sec. IVB comes into play. At higher 
temperatures the tunneling band is activated, so that 
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S w = InL + f(j ! ^) ; m 1 ' 23 L~ 2 « T « H/S . (76) 

Here In L + F(x) is the entropy of a one dimensional tight binding band of L sites. By noting 
the form of the density of states one can see that for small x, F(x) ~ (1/2) lnx + Const and 
as x — > oo, F(x) — > 0. The upper limit of this regime is defined so that essentially only the 
lowest subband of Sec. IVB is activated. Next, there is a regime when the temperature is 
high enough that all states of these subbands are equally populated, but only one domain 
wall exists. In this discussion it is assumed that H ^ J/S, so that the band of particle 
states is well separated from the states with more than one wall. In this regime one has 

S w = ln(2SX) +F(jj) , H/S^T < J/S . (77) 

In this case, F(x) — > for x — > oo and F(x) ~ \n(T/2H) + 1 for x — > 0. Finally, we have 
the classical regime when T ^> J/S. In that regime any effects due to the fact that S x (i) 
and S z (i) do not commute with one another becomes unimportant. In that regime then, 
the partition function for the quantum model of Eq. (^) is the same [up to corrections of 
relative order J/ (ST)] as that of the classical model in which the operators are interpreted 
as classical variables according to Eq. (|T6"D. 

The effect of a random potential is somewhat different. In the present context a random 
potential would be created by a field which independently for each site assumes random 
values from a distribution of width Vr. It is knownil that in the presence of a random 
potential all states are localized. That means that in the presence of a random field the 
eigenfunctions are localized and therefore the width of the wall in the ground state remains 
finite in the limit L — > oo. Nevertheless, the width of the wall, W, can be expressed in terms 
of scaling functions of the same variables as describe the thermodynamics except that here 
Vr replaces T. Thus 

W ~ L ; Vr < H2 1 - 2S L~ 2 , (78) 

and 
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W / a = F (jj2^s) ; H2 1 ~ 2S L~ 2 « V R « H/S . (79) 
Since we know that W can not be less than the intrinsic width, we see that F(oo) = 2. 
Similarly, we have 

W/a = F(jf); H/S « V R « J/S . (80) 

In this case -F(O) = 2 and F(oo) = 1, the intrinsic wall width in this regime. For any larger 
Vr, W/a remains equal to unity. 

VI. SUMMARY 

We may summarize our conclusions for the Ising model in a transverse field as follows. 

1. Although an exact analysis of domain walls in the transverse Ising model can in 
principle be accomplished using the Jordan- Wigner transformation, it is much 
simpler to consider the nature of degenerate perturbation theory in the presence 
of "up-down" boundary conditions. Within degenerate perturbation theory for 
S = 1/2 and S = 1 the Hamiltonian governing the position of the center of the 
domain wall is isomorphic with a hopping model with a hopping matrix element 
t ~ H. We verified this in detail with analysis of the scaling of the energy levels 
with chain length, L. 

2. For spins |, 1 and | the ground state profiles when spins at opposite ends of a 
chain of length L are fixed to be antiparallel show relatively little dependence on 
the transverse field h (as long as h is less than the critical value, h c , above which 
ordering in S z is destroyed). For small h the profiles at large but finite S are 
not very different from those at small S. In all cases, the profiles are quite well 
obtained from the hopping model. 

3. The profiles in the oscillatory phase show much larger finite size corrections than 
those in the ferromagnetic phase. In the limit L — > oo the profiles in the oscilla- 
tory phase may still show some oscillatory behavior, unlike in the ferromagnetic 
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phase. Also, our results for the profiles show an interesting finite-size effect in 
which spatial parity appears to be broken. 

4. In the spectrum of the oscillatory phase at h = one finds an unusual pattern 
for the gaps between the ground state and the first few excited states. While for 
L even, the energies form doublets, each with a splitting which scales with L as 
L~ 3 , and with gaps between the doublets which scale as L~ 2 . For L odd, the 
energies are doubly degenerate and the gaps scale with L as L~ 2 . This behavior 
was explained in terms of dimer excitations which occur in the limit h — > when 

5. For large S the quantum transverse Ising model the Hamiltonian describing the 
position of the center of the domain wall gives rise to band splittings analogous 
to those found for an electron in a periodic potential^. The lowest sub-band is 
described by a hopping model in which the center of the domain wall tunnels 
from one site to a neighboring sites. The associated tunneling matrix element is 
of order hS exp(-aS), where a = 2 In 2. It would be interesting to recover this 
result within a field theory for one space and one time dimension. 

6. For the classical (S = oo) model the domain wall is extremely narrow for small 
h. We have determined the domain wall energy for the case when i) the center of 
the domain wall is at a lattice site and ii) the center of the domain wall is midway 
between two lattice sites. The difference between these energies is the pinning 
energy which must be exceeded to move the domain wall through the lattice. For 
the quantum model we find this pinning energy vanishes as L — > oo due to the 
fact that the width of the domain wall diverges in this limit. 

7. In attempting to relate quantum problems to classical problems it is necessary to 
define an intrinsic width, Wi, of a quantum domain wall. We have defined Wi as 
the minimum possible width attainable using wave functions from a manifold M.. 
At temperatures large compared to the band but small compared to the activation 
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energy for creation of an additional wall, the manifold M. is the hopping band. 
In that case, the intrinsic width of the quantum system is essentially equal to the 
width of the analogous classical domain wall. 
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TABLES 

TABLE I. Finite-size estimates of the exponent 6 for the spin-1 model for two values of h. 
The line labeled oo gives the L — > oo estimates and the numbers in brackets give the estimated 
uncertainty in the last given digit. 
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TABLE II. Extrapolated estimates for the finite-size scaling amplitudes a(i) for S = ^ i n the 
ferromagnetic phase F. The last column gives the amplitude ratios r(i) = a(i)/a(l) from Eq. (67). 
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TABLE III. Extrapolated estimates for the finite-size amplitudes a{i) for spin 5=1. 
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TABLE IV. Estimates for lO 5 ^ and lO 5 ^, giving upper bounds for the wall energy differences 
2\B-B'\ at ft = 0.1 
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TABLE V. Finite-size data for the exponent 9 obtained for S = \ with h = and e = 0.02 
and the extrapolated limit for L — > oo. Estimates obtained for L even and odd, respectively, are 
shown separately. 
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TABLE VI. Extrapolated finite-size amplitudes a(i), i = 1, 2, . . . , 6 for S = \ with h = and 



for several values of e. For a given e, the first (second) line corresponds to the estimates found for 
L even (odd). 



e 


i 


z 


Q 
O 


A 

4 


K 



D 


0.02 





1.2033(20) 


1.22(4) 


3.191(20) 


3.243(15) 


5.97(10) 







1.203(3) 


1.203(3) 


3.191(15) 


3.191(15) 


5.975(20) 


0.10 





6.553(15) 


6.58(3) 


17.26(15) 


17.59(10) 


32.3(3) 







6.555(20) 


6.555(20) 


17.32(8) 


17.32(8) 


32.48(10) 


0.50 





57.2(9) 


58.8(6) 


153(4) 


158(3) 


311(9) 







58.3(3) 
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FIGURE CAPTIONS 



Fig. 1 Phase diagram of the "yz" model of Eq. ([|) for spin S = |. The disordered (D), 
ferromagnetic (F) and oscillatory (0) phases are indicated. 

Fig. 2 Local magnetization profile m(r) for the spin-| model with h = 0.1 and e = for 
system sizes L — 2 ... 16. 

Fig. 3 Comparison of the magnetization profiles m(r) for the spin-| (boxes), the spin-1 (full 
circles) and the spin-| (open circles) models. The full curve gives the profile of Eq. flllf) 
which is correct to first order in h. The data are for e = and A) h = 0.1 and B) 
h = 0.5. 

Fig. 4 Magnetization profile m(r) for three values of the transverse field h and for e = and 

s = i 

Fig. 5 States in the ground manifold M.. Top: the states \n) introduced for spin 1/2. Bottom: 
the states |n) introduced for the case S ^> 1. In this case the nth spin has S x = S. 

Fig. 6 Magnetization profile of the oscillatory phase of the spin-| model with e = 0.5 and 
h = 0.01, for A) L even and B) L odd. 

Fig. 7 Comparison of the magnetization profiles of the spin-| model with h — > and the 
values of e indicated (data points) and with h = h e{i (e) = yfe and e = (curves), for 
L = 15. 

Fig. 8 Density of states (per site) D(E) as a function of energy E (relative to the center of 
the manifold) within the manifold M. for 5* = 3. The middle subband contains 2L + 1 
states. The other subbands contain L states. The asymmetry of the subbands rapidly 
decreases as S increases. Inside each subband, near its edge at Ei, D(E) ~ \E — E^ 1 / 2 
for L — > oo. 
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Fig. 9 Schematic structure of the low- lying states for spin S — \, h — and e 7^ 0, for A) L 
even and B) L odd. 

Fig. 10 Spectrum of H yz for spin S — \ and e = 0.25 as a function of h for L = 10 sites. The 
two vertical lines indicate the location of the transitions between the O /F phases and 
the F/D phases, respectively, for L — > 00. 
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